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Abstract 

Geometric Singular Perturbation Theory (GSPT) and Conley Index Theory are two 
powerful techniques to analyze dynamical systems. Conley already realized that using 
his index is easier for singular perturbation problems. In this paper, we will revisit 
Conley's results and prove that the GSPT technique of Fenichel Normal Form can 
be used to simplify the application of Conley index techniques even further. We also 
hope that our results provide a better bridge between the different fields. Furthermore 
we show how to interpret Conley's conditions in terms of averaging. The result are 
illustrated by the two-dimensional van der Pol equation and by a three-dimensional 
Morris-Lecar model. 
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1 Introduction 

We start by outlining the context and results of the paper. Singular perturbation theory 
often involves a distingushed parameter, usually denoted e, which is assumed to be small 
and positive. For e = 0, the dynamical system is "degenerate" but it might also be easier 
to analyze. Ordinary differential equations (ODEs) with two different time scales, so-called 
fast-slow systems, are a very important class of singular perturbation problems. In this 
context, the parameter e describes the separation of time scales. One successful strategy to 
analyze fast-slow systems is to understand the case e = and using this knowledge to try 
to prove perturbation results for e > sufficiently small; many geometric and asymptotic 
methods follow this pattern. 



Conley Index Theory can be applied to wide classes of dynamical systems which do 
not have to be singularly perturbed. The goal is to convert the problem into an alge- 
braic/topological question. Conley already proved that this approach can be substanially 
simplified for singular perturbation problems. Mischaikow and co-workers developed Con- 
ley's index theory for fast-slow systems recently even further. The problem with the broad 
applicability of the theory is that it is still very technical. It has been demonstrated in 
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low- dimensional examples that easily applicable geometric versions of the theory should be 
possible. One possible generalization to higher dimensions has been proposed but requires 
a rather complicated and lengthy topological construction. 

The goal of the current paper is to start building a simpler bridge between GSPT and 
Conley Index Theory; in particular, we are going to show that a fundamental result due to 
Conley can be simplified using Fenichel Normal Form. The result we prove shows that if we 
isolate an invariant set by a suitable neighbourhood N, then N is of the form required by 
Conley Index Theory if the slow motion for e = is transverse to the boundary of N. We 
also investigate the case when the fast motion of the system has periodic orbits and discuss 
further generalizations. We hope that the new results will be of interest from the viewpoint 
of GSPT as well as Conley Index Theory. 

The paper is structured as follows. In Section 2 we describe the background from fast- 
slow systems theory and in Section 3 a similar exposition is given for Conley Index Theory 
with a focus on the application to fast-slow systems. Both introductions are a little more 
detailed than strictly necessary to accomodate the two different perspectives with respect to 
background knowledge. In Section 4 we prove the main result for equilibrium points of the 
fast motion. In Section 5 a result for the case of periodic orbits is given and a generalization 
to any bounded invariant set is outlined. Furthermore the role of averaging is explained. In 
Section 6 we demonstrate the applicability of the result by two examples. 



2 Fast- Slow Systems 

Several viewpoints have influenced the development of multiple time scale or fast-slow sys- 
tems starting with asymptotic analysis [261 IB] using techniques like matched asymptotic 
expansions [20, 23J. A geometric theory focusing on invariant manifolds was developed 
[HI HS1 E2] which is now commonly known as Fenichel theory due to Fenichel's seminal work 
[9]. There was also significant influence by a group using nonstandard analysis 0Q]. 

We shall focus on the geometric viewpoint in this paper. The term "Geometric Singular 
Perturbation Theory" (GSPT) is used to encompass Fenichel theory and further geometric 
methods developed over the last three decades in the context of multiple time scale problems. 
The general formulation of a fast-slow system of ordinary differential equations (ODEs) is 

ex = ef = f(x,y,e), m 

dy i \ \ l ) 

y = ■£ = g{x,y,e), 

where (x, y) G R m x R n and e is a small parameter < e < 1 representing the ratio of time 
scales. The functions / : R m x R n x R -»■ R m and g : R m x R" ■ x R R" will be assumed 
to be sufficiently smooth. The variables x are fast and the variables y are slow and we can 
change in ([1]) from the slow time scale r to the fast time scale t = r/e which yields: 

x ' = % = f(x,y,e), (] 
y> = i = eg(x,y,e). {A) 
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We will also denote the vector field by z' = F(z) where F = (/, eg) and z = (x, y). The 
first major idea to analyze (JTJ-flSJ) is to consider the singular limit as e — > 0. 

Definition 2.1. Setting e = in (J2J gives 

y' = o, ^ 

which is system of ODEs parametrized by the slow variables y. We call (|3J) the /as£ subsystem 
or /ayer equations. The associated flow is called the /as£ flow. 

Definition 2.2. Considering the singular limit e = for (JTJ yields: 

= f(x,y,0), ,^ 

y = g(x,y,o). 

System (j4]) is a differential-algebraic equation (DAE) called slow subsystem or reduced system. 
The associated flow is called the slow flow. 

One goal of GSPT is to use the fast and slow subsystems to understand the dynamics of 
the full system ©-© for e > 0. 

Definition 2.3. The algebraic constraint of fll]) defines the critical manifold 

C := {(x,y) eR m xR n \f(x,y,0) = 0}. 

Note that it is possible that C is not an actual manifold [22] but we shall not consider 
this case here. The points in C are equilibrium points for the fast subsystem (jHJ). 

Example 2.4. Consider the following very simple planar fast-slow system 

ex = y-x 2 , f . 
y = -1. W 

The critical manifold C = {(x, y) G M? : y = x 2 } is a parabola. Observe that the slow flow 
on C is y = — 1 so that under this flow any initial condition on C will "flow down" to the 
origin (x,y) = (0,0); see single arrows in Figure Ufa). The fast subsystem is x' = y — x 2 
which has one stable equilibrium and one unstable one for y > 0, a saddle- node (or fold) 
bifurcation for y = and no equilibria for y < 0; the flow is indicated by double arrows in 
Figure QJa). 



Definition 2.5. A subset McCis called normally hyperbolic if the mxm matrix (D x f)(p) 
of first partial derivatives with respect to the fast variables has no eigenvalues with zero real 
part for all p G M; this condition is equivalent to requiring that points p G S are hyperbolic 
equilibria of the fast subsystem (EJ). 
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Figure 1: (a) Critical manifold C (dark grey) of (JSJ); the fast and slow flows are indicated 
by double and single arrows respectively, (b) C (dark grey) and a slow manifold M e (light 
grey) obtained from Fenichel's Theorem 12.71 are shown together with a trajectory 7 of (jSJ) 
starting at 7(0) = (1,0.4) at parameter value e = 0.05. Observe that 7 quickly approaches 
M e , then tracks it (actually 0(e~ K ^ e )-c\ose) until the fast-slow structure breaks down near 
the fold point (x, y) = (0,0). 



We call a normally hyberbolic subset M attracting if all eigenvalues of (D x f)(p) have 
negative real parts for p e M; similarly M is called repelling if all eigenvalues have positive 
real parts. If M is normally hyperbolic and neither attracting nor repelling we say it is 
of saddle-type. A typical class of points where normal hyperbolicity fails are fold points. 
They are defined as the points where the critical manifold C is locally parabolic with respect 
to the fast directions. In other words, at a fold point p* one requires that /(p*,0) = 
and that (D x f)(p*,Q) is of rank m — 1 with left and right null vectors w and v, such that 
w ■ [{D xx f)(p){v,v)\ ± and w ■ [{D y f){p)} ± 0. 

Example 2.6. (Example 12.41 continued) The critical manifold C = {y = x 2 } splits into one 
repelling part, one attracting part and a fold point at the origin: 

C = C z U{(0,0)}UC r 

where C\ = C fl {x < 0} is repelling, C r = C PI {x > 0} is attracting and (x, y) = (0, 0) is a 
fold point which is easily verified since 

df d 2 f 
DJ = jf = -2x and D xx f = —L = -2. 
ox ox z 

We continue with the general case. If (D x f) has maximal rank the implicit function 
theorem applied to f(x,y,0) = locally provides a function h(y) = x so that C can be 
expressed as a graph. Hence the slow subsystem (TJj) can be more succinctly expressed as: 

y = g(h(y),y,o) (6) 

We shall also refer to the flow induced by ([6]) as slow flow. To relate the dynamics of the 
slow flow to the dynamics of the full system for e > the next theorem is of fundamental 
importance. 
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Theorem 2.7 (Fenichel's Theorem, [9, 18J). Suppose M = M is a compact normally hyper- 
bolic submanifold of the critical manifold C . Then for e > sufficiently small the following 
holds: 

(Fl) There exists a locally invariant manifold M e diffeomorphic to Mq. Local invariance 
means that M E can have boundaries through which trajectories enter or leave. 

(F2) M t has a distance 0(e) from M Q . 

(F3) The flow on M f converges to the slow flow as e — > 0. 

(F4) M e is C r -smooth for any r < oo (as long as f,g G C°°). 

(F5) M t is normally hyperbolic and has the same stability properties with respect to the fast 
variables as M (attracting, repelling or saddle-type). 

(F6) For fixed e > 0, M e is usually not unique but all manifolds satisfying (F1)-(F5) lie at 
a Hausdorff distance 0(e K l e ) from each other for some K > 0, K = 0(1). 

We call a manifold M t a slow manifold. Note that all asymptotic notation refers to e — > 0. 
The same conclusions as for Mq hold (locally) for its stable and unstable manifolds: 

W S (M ) = |J W s (p), W U (M ) = |J W u (p) 
peM peM 

where we view points p G Mq as equilibria of the fast subsystem. 

Figure [H(b) shows an typical scenario where Fenichel's Theorem applies; there we picked 
a compact submanifold Mq C C r and obtained an associated slow manifold. In addition 
to Fenichel's Theorem we can also find coordinate changes the simplify a fast-slow system 
considerably near a critical manifold. 

Theorem 2.8 [Fenichel Normal Form, [HI EH])- Suppose the origin G C is a normally 
hyberbolic point with m u unstable and m s stable fast directions; choose a sufficiently small 
compact normally hperbolic subset Mq C C containing the origin. Then there exists a smooth 
invertible coordinate change (x, y) i— >■ (a, b, v) G M m " x R ms x M n in a neighbourhood of so 
that a fast-slow system (j5J) can be written as 

a' = A(a,b,v,e)a, 

b' = T(a,b,v,e)b, (7) 
v' = e(m(v,e) + H(a,b,v,e)ab), 

where A, T are matrix-valued functions. A has m u positive and Y has m u negative eigenval- 
ues. H is bilinear and given in coordinates by 

Hi(a,b,v,e)ab = H irs a r b s . (8) 

r=l s=l 
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Figure 2: Illustration of Theorem 12.81 



The situation is illustrated in Figure EJ The manifold Mq perturbs to a slow manifold M e 
by Fenichel's Theorem. Then this slow manifold is "straightened" together with its stable 
and unstable manifolds that become coordinate planes. 

As a major overall conlusion of Fenichel Theory we get that the flow near a normally 
hyperbolic critical manifold is "completely determined" by the singular limit systems for 
e = 0. We will show in Section H] how this result reappears for the Conley index theory of 
fast-slow systems. 

3 Conley Index Theory 

In this section we describe the basic constructions from Conley Index Theory and how these 
have been adapted to fast-slow systems. When Conley [3] studied these techniques that now 
bear his name it seems possible [TTJ that he also had applications to singular perturbation 
problems in mind. The idea of the theory is to convert a dynamical problem (e.g. "Does 
this dynamical system have a heteroclinic orbit") into an algebraic problem (e.g. "What 
is the structure of a matrix?"). Several successful applications exist; see [23] for a recent 
survey. We are going to outline only the basic techniques of the theory focusing on isolating 
neighbourhoods for fast-slow ODEs. In this context, the current theory can be found in 
[3j |25j [TTJ [TOj [2TJ . For a more detailed introduction in the case of general finite-dimensional 
dynamical systems we refer to (23J; the infinite-dimensional case is considered in [29] . 

Let : 1R x M fe — > M fc be a flow with <p = <f)(t,z). A compact set N C R k is called an 
isolating neighbourhood if 

Inv(JV, 0) := {z e M fc |0(M, z) C iV} C int(JV) 

where int(iV) denotes the interior of N . If we set S := Inv(iV, <ft) then S is called an isolated 
invariant set. The situation is illustrated in Figure [3] where Inv(iV, 0) is an unstable node. 
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Figure 3: An isolating neighbourhood N (shaded disk) of an unstable node q is shown. The 
node q is the invariant set of of N i.e. Inv(iV, <p) = {q}. The boundary of N (dashed circle) 
is denoted by L and (N, L) is easily seen to form an index pair. 



Definition 3.1. Let S be an isolated invariant set. A pair of compact sets (N,L) with 
L C iV is called an index pair for S if the following conditions hold: 

(a) S = Inv(cl(iV — L)) and N — L is a neighbourhood of S. 

(b) L is positively invariant in N i.e. for any z G L and 0([O, t],z) G N then 0([O, t],z) C L. 

(c) L is an exit set for iV i.e. for any z G N and £i > such that <p(ti, z) G" N then there 
exists to G [0,£i] for which </>([0, t ], z) C N and 4>(t ,z) G L. 

We define the Conley index of S as 



where if* is relative homology [TBI [3D] . Note that an alternative way to define the Conley 
index would be to consider cohomology and set CH*(S) := H*(N,L). The strategy to use 
the Conley index for dynamical systems usually proceeds along the following lines: 

(51) Find an isolating neighbourhood N. 

(52) Determine an index pair (N,L). 

(53) Calculate the Conley index. 

(54) Use the calculation to prove a result about Inv(iV, 0). 

We will focus on (S1)-(S2) in the context of fast-slow systems. The main question is whether 
we can use the fast-slow structure to find an index pair. Re-writing a general fast-slow 
system (|2J) on the fast time scale with z — (x,y) will be convient 



We denote the flow of © by <p e : R x R m+n -> R m+n . Observe that O is the flow of the fast 
subsystem ([3]). It is easy to see that if N is an isolating neighbourhood for 0o then it is also 
an isolating neighbourhood for <p t . The problem is that usually iV will not be an isolating 
neighbourhood for e = but it still can be an isolating neighbourhood e > 0. 



CH*{S) :=H,(N,L) 



j 




(9) 



i=i 
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Definition 3.2. A compact set N C M m+ra is called a singular isolating neighbourhood if 
N is not an isolating neighbourhood for <p but there exists e such that N is an isolating 
neighbourhood for <p € with e G (0, e]. 

The next example illustrates, without proof, a singular isolating neighbourhood in a 
fast-slow system. 

Example 3.3. A time reversed version of Van der Pol's [5J E] equation is 

x x /3 x y, (ic\\ 
y' = ex. 

We shall consider the Van der Pol equation in more detail in Section 16.11 For now it is 
useful to look ahead to Figure [6] that shows the critical manifold Co of f JTOj) and an orbit for 
e = composed of fast and slow subsystem trajectories. To prove that this orbit perturbs we 
want to construct a singular isolating neighbourhood; the dashed lines in Figure [6] indicate 
a possible guess for a such a neighbourhood N. In Section IBTTl we are going to prove that N 
is a singular isolating neighbourhood. For now it is important to observe that N is not an 
isolating neighbourhood for e = on the fast time scale (check it!). 

To check whether a compact set is a singular isolating neighbourhood we temporalily 
decide to define the complications away. 

Definition 3.4. Let iV be a compact set and let z G Inv(iV, 0o) =: S. We say that z is a 
slow exit (entrance) point if there is a neighbourhood U of z and an e > such that for all 
e G (0, e] there is a time T(e, U) > (T(e, U) < 0) such that 

<P t (T(e,U),U)nU = ®. 

Let S~ (S + ) denote the set of slow exit (entrance) points. Furthermore define the following 

sets 

S d :=SndN and S^:=S d nS ± . 

To understand what slow exit and entrance points are and what characterizes them is 
the main goal of this paper and should be clear after Section HJ Obviously Definition 13.41 
'cheats' by prescribing the dynamics of the slow motion under perturbation. Therefore the 
next result is very easy to prove. 

Theorem 3.5. If Inv(N, O ) H dN consists of slow exit and entrance points then N is a 
singular isolating neighbourhood i.e. it is an isolating neighbourhood for the full fast-slow 
system for sufficiently small e > 0. 

Hence we have reduced the problem to characterizing slow exit /entrance points in a more 
computable way. The next definition provides a technical notion which will be necessary for 
this task. 

Definition 3.6. The average of a function h on S C M. m+n , denoted Avg(h, S), is the limit 
as T — > oo of the set of numbers 

Jo h ^°( s,z ^ ds z e s 

We say that h has strictly positive averages on S if Avg(/i, S) C (0, oo). 
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In a seminal paper [3] Conley was able to give computable conditions for slow exit and 
entry points. 

Theorem 3.7. (fE[) A point z G S is a slow exit point if there exists a compact set K z C S 
invariant under <pQ, a neighbourhood U z of the chain recurrent set 1Z(K Z ) of K z , an e > 
and a function I : cl(U z ) x [0, e] — > R such that the following conditions are satisfied: 

(a) u}(z,4>o) = HieR c U{^o(s, z) : s > t}) C K z ; here cl{.) denotes closure. 

(b) I is of the form l(w, e) = lo(w) + eli(w) + . . . + eHj(w). 

(c) If Lq = {w\lo(w) = 0} then K z Dcl(U z ) = SnL nd(U z ) and furthermore lo\snd(u x ) < 0. 

(d) Let Gj(w) = V 2 Z H • Fj(w) + V z h(w) ■ i^--i(w) + . . . + V z ?» • F (w). Then for 
some k, Gk = if k < j and Gj has strictly positive averages on 1Z{K Z ). 

A point is a slow entrance point if the same conditions hold under reversal of time. Points 
that satisfy (a)-(d) (satisfy the conditions under time reversal) are called C-slow exit (en- 
trance) points. The compact set K z is called a slow exit guide. 

We note that the function / should be viewed as a Lyapunov-type function for the dy- 
namics near the slow exit /entrance point. In [25] the authors claim that "the only dynamics 
which plays a role in the calculations (for Theorem 13.71) is that of 0o" • Formally this is not 
the case since higher-order terms Fj for j > do enter crucially in (d). But as we shall see 
in the next section, the idea was intuitively correct. In fact, one should state that only the 
fast and slow flows of the singular subsystems play a role in the calculations. 

Using Theorem 13.71 we can often identify an isolating neighbourhood for a fast-slow 
system. Mischaikow, Mrozek and Reineck [25] give an analogous construction for index 
pairs. 

Definition 3.8. A pair of compact sets (N, L) with L C iV is called a singular index pair 
if cl(N — L) is a singular isolating neighbourhood and there exists an e > such that for all 
ee(0,e] 

H*(N, L) = C#*(Inv(cl(iV - L), <f> e )). 

The singular index pair should be characterized by similar conditions as the usual index 
pair described in Definition l3.il From the exit set requirement we know that L has to contain 
the immediate exit set of N 

N~ := {z G dN\<p ((0,t) 7 z) £ N for all t > 0}. 

Regarding positive invariance, it turns out that give Y <Z N one has to consider the push- 
forward set in iV under the flow Y defined by 

p(Y,N,(f) ) := {z G N\3w G Y,t > s.t. O ([O, t], w) C N,(j) (t,w) = z). 
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Basically p(Y, N, O ) consists of points in N that can be reached from Y by a positive 
trajectory in N; observe that by construction we must have Y C p(Y,N,<p ). In addition, 
we also must consider a special version of the unstable manifold of a point lying in N 

W%(Y) := {z G iV|0o((-oo,O),z) C N and a{z,<j> Q ) C Y}. 

Again we observe that Y C W^(Y). Before we can state the theorem about characterizing 
singular index pairs, one last definition is needed. 

Definition 3.9. A slow entrance point z is called a strict slow entrance point if there exists 
a neighbourhood V of z and an e > such that if v G V fl iV and e G (0, e] then there exists 
a time t v (e) such that 

(j) e ([0,t v (e)],v) C N. 
The set of strict slow entrance points will be denoted by Sq + . 

Theorem 3.10. (l2h^ ) Let N be a singular isolating neighbourhood. Assume 

(A) Sq consists of C-slow exit points. 

(B) S d cS+ + US d . 

(C) - S 9 ) n cl(N-) = 0. 
For each z G Sq, let K z be a slow exit guide for z. Define 

L := p(cl(N~), N, O ) U W£ | |J H{K Z ) 

\zeSg 

If L is closed then (N,L) is a singular index pair. 

Observe that for Theorem 13.71 and Theorem 13.101 it is crucial to determine which points 
are slow exit /entrance points. The conditions (a)-(d) given in Theorem 13.71 are complicated. 
The goal of this paper is to simplify these conditions. 

4 Equilibrium Exit Points 

Let iV be a compact set and let z = (xq, yo) G Inv(iV, <p ) =: S where 0o denotes the flow of 
the fast subsystem. Let C denote the critical manifold. We make the following assumptions: 

(Al) z G c n ON. 

(A2) C is a normally hyperbolic manifold at z$ and locally given as a graph x = h(y). 
(A3) dN is smooth and parallel to the fast fibers near zq. 

(A4) The slow flow y = g(h(y),y) is transverse to dN near z . Let n denote the outward 
unit normal to iV at Zq; there are two cases: 




10 



(A4.1) n ■ (0, g(h(y ), y )) > 0, slow flow directed outward near z . 
(A4.2) n ■ (0, g(h(yo),yo)) < 0, slow flow directed inward near z . 

Theorem 4.1. Under conditions (A1)-(A4) the point zq is a slow exit/entrance point. If 
(A4-1) holds we have a slow exit point and for (AJ^.2) we get a slow entrance point. 

Remark: Theorem 14.11 only requires knowledge about the fast and slow subsystems for 
e = 0. This justifies more clearly than Theorem 13.71 that a special Conley index theory for 
fast-slow systems is possible. 

Proof. Suppose without loss of generality that (A4.1) holds so that we are trying to show 
that Zq G C is a slow exit point. We also work in a sufficiently small neighbourhood of z 
for the rest of the proof. The goal is to verify the conditions (a)-(d) of Theorem 13.71 Using 
(A2) we apply Fenichel's Normal Form Theorem 12.81 to re-write the fast-slow system near 
as 

x' = n(x,y,e)x, , n - 
y' = t(m(y,e) + H(x,y,e)x), 

where H is bilinear as described in (JS} and we have re-written the Fenichel coordinates as 
x = (a, b) e W 1 and v = y e W 1 with Q = (A, T). Due to (A3), we obtain that dN is locally 
given by {y = 0} and also locally we have N = {yi < for alH = 1, . . . , n}; see Figure HI 




Figure 4: Illustration of the situation near a slow exit point at the origin in IR 3 . The compact 
set N is locally given by {y < 0}. The outer normal vector e± to dN is also shown. The 
slow flow will point along this normal vector. The parabolic surface Lq is the zero set 
l(x,y) = = l (x,y) . 



Note that (A4.1) implies that near we can rectify the slow flow so that ffTTj) becomes: 

x' = Sl(x,y,e)x, , l2 ^ 
y' = e(ei + H(x,y,e)x), 
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where e\ = (1, 0, . . . , 0) G M. n . As a slow exit guide set K = K := {0} and observe that 
since G C it is an equilibrium point for the fast subsystem. Therefore u(0, 0o) — {0} and 
so u>(0, 0) = K which verifies (a). Define the function I by 

in 

i( x >y) = Vi- ^2( x j) 2 = kfay). 

Notice that I = l and so we find that (see also Figure 0J 

L = {(x,y)\l (x,y) = 0} = = ^(a^) 2 j . 

Let U be a sufficiently small neighbourhood around then U H K — {0}. We also have 
locally S n c\(U) = N n C and therefore 

5" n c\(u) n Lq = {o} 

since = Y^j=i(. x i) 2 holds if and only if Xj = for all j = 1,2, ... , m. Obviously lo\snc\(u) < 
and so (b)-(c) hold. For the last step, observe that 7^.(0) = {0} and hence we have to verify 
that condition (d) holds at the origin i.e. 

Gj(z) = VM*) ■ Fj{z) + V z h{z) ■ Fj^(z) + ... + V z lj{z) ■ F {z) 

satisfies that for some k, Gfe(0) = if k < j and Gj(0) > 0. We compute the gradient of lo 

VJ (x, y) = (-2a;i, -2x2, • • • , -2x m , 1, 0, . . . , 0). 

Since 0o describes the fast flow, the first term F in for the normal form (fT2"|) is given by 

F (x,y) = (n(x,y,0)x,0) T . 

This gives that V 2 /o(0) ■ -F (0) = ■ Q(x, y, 0)x + e% ■ = 0. Hence Go(0) is identically zero. 
Next, we show that Gi(0) is positive. We have 

Gt{z) = VMz) ■ F i(z) = (-2^1, -2x 2 , . . . , -2x m , 1, 0, . . . , 0) • (0, e x + H(x, y, e)x) T . 

Since H(0, 0, e)0 = we immediately get Gi(0) = e\ ■ [t\) T = 1 > verifying (d). Therefore 
the original point zq is a slow exit point. □ 

The only condition that does not seem not quite natural for Theorem 14.11 is (A3). To 
illustrate that it is necessary consider the following example. 

Example 4.2. Consider a fast-slow system with (x,y) G M 2 given by 

x' = x, 
y = e. 

The solution is given by (x(t),y(t)) = (x(0)e*, y(0) + et). Fix some m > and let iV = 
{{x,y) G M?\y < mx} locally near i.e. we truncate iV outside a suitable neighbourhood to 
make it compact. Now the origin is not a slow exit point although (Al)-(A2) and (A4.1) hold. 
Indeed, pick a neighbourhood U of then there is (x(0), y(O)) G U such that x(0),y(0) > 0. 
For e > sufficiently small we can easily assure that y(0) + et < mx(0)e* for all t > such 
that the trajectory starting at (x(0),y(0)) stays in N. 
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Example 14.21 also indicates that the problem should not occur for attracting critical/slow 
manifolds. 

Proposition 4.3. Suppose (A1),(A2) and (A4) hold. Furthermore assume that dN is locally 
linear and has an angle of order 0(1) to C at zq and that C is attracting at zq. Then zq is 
a slow exit/ entrance point. If (A4-1) holds we have a slow exit point and for (A4-2) we get 
a slow entrance point. 

Proof. In this case the proof is much simpler and we do not need Conley's Theorem 13.71 
Again we can restrict without loss of generality to the case (A4.1). We will work under the 
assumption for the rest of this proof that e > has been chosen so that Fenichel Theory 
applies. Applying Fenichel' s Normal Form Theorem as in Theorem 14. II gives 

x' = Q(x,y,e)x, , 13 . 
y' = e( i m(y,e)+H( l x,y,e)x), 

where now Q(0, 0, 0) has m negative eigenvalues. Let U be a small neighbourhood around the 
origin. By Fenichel's Theorem 12 . 71 ( F 1 ) there exists a slow manifold C e . Let 7 be a trajectory 
with an initial condition in U. Since Q(0, 0, 0) has m negative eigenvalues Fenichel's Theorem 
(F5) shows that C t is attracting. Therefore 7 gets attracted exponentially to C t or lies in 
C e . Observe that dN is not tangent to C e as it has an 0(1) angle to Cq. By (A4.1) and 
Fenichel's Theorem (F3) we find that 7 must leave N after a time T 7 (e). Since U is bounded 
we can take the maximum of all times over cl(U) 

T(e,U) := max TJe). 

7(0)ecl(J7) 1 

This verifies Definition 13.41 □ 

Unfortunately Proposition 14.31 is rarely helpful. One reason is that often it is convenient 
to make the critical manifold repelling near slow exit points; see examples in Section The 
main reason is that in many important cases critical manifolds of saddle- type appear [13]. 
In fact, one of the most well-known examples, the 3D FitzHugh-Nagumo equation, has two 
fast variables and one slow variable with a critical/slow manifold of saddle type [TJ1 [15]. 

5 Periodic Orbit Exit Points 

In this section we shall not aim for the most general results but show some characterizations 
of slow exit points in the case of periodic orbits. We restrict to the case of fast-slow systems 
in IR 3 with two fast variables i.e. (x, y) G M 2 x M and 

z'i = fi(x,y), 

x' 2 = f 2 (x,y), (14) 
V = eg(x,y). 

Let 7y(t) G M 2 denote a periodic orbit for the fast subsystem with period T y so that 

7,(0) = ly (T y ), 1 ' y (t) = f( ly (t),y). 

Let iV := [-K, K] 2 x [-K, 0] C M 3 for K > so that the following assumptions hold (see 
Figure ED : 
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(Bl) There exists family of hyperbolic periodic orbits { 7y } for y G [—80, 8 ] for some 8 > 
in the fast subsystem. 

(B2) {l y }ye[-8 , 0] C N and 7 o C mt([-K,K] 2 x {0}). 




N 



Figure 5: Sketch of the situation near periodic orbit 70 C dN in the fast subsystem. The 
parabolic surface L given by the zero set l (x,y) = is defined by rotating the given 
parabola along 70. 



Proposition 5.1. Suppose (B1)-(B2) hold and assume that 

i I ° g( lo (s),0)ds>0 (15) 
J o Jo 

where g is as given in ( fT4"l) . Then all points in 7 o are slow exit points for N. 

Proof. The proof is very similar to the argument for Theorem 14. 11 Let z e = (x e , 0) G 7 o 
be any point in the periodic orbit contained in dN. Observe that u(x e , 0o) — 7o an d let 
K = 7o = 1Z(K). Let U be an annular neighbourhood of K contained in N, for example we 
can set 



U = <i (x, y) G IR 3 : min||z 7 — (s,y) T || 2 < 5\ for G 7y with y G [—82,82] 



for 81,82 > sufficiently small, let 7r : IR 3 — > M 3 denote the orthogonal projection onto 7q. 
Now define 

2 

l(z) = l(x,y) = l (x,y) :=y- ^(xj - M z )j) 2 

where the subscript j indicates the Xj-coordinate of a point. In the notation of Theorem 
13.71 we easily check that (c) holds and we also observe that for (d) we have G (z) = since 
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V z l = (0,0, 1) T and F = (/i,/ 2 ,0). We also find that G x = (0,0, g(z)) ■ (0,0, 1) T and on 
K we indeed have 

r-To 



^ / ° g(lo(s),0)ds>0 
I o Jo 



which verifies (d) and shows that z e is a slow exit point. Noting that z e was arbitrary on 70 
finishes the proof. □ 

Note that Proposition 15.11 has the rather obvious interpretation that a point is a slow 
exit point for a periodic orbit of the fast subsystem if it lies on the boundary of the compact 
set N and the average slow drift moves it outside of N. It is more interesting to re-interpret 
the condition (fT5l) . 

We want to deal with families of periodic orbits in the fast subsystem. For normally 
hyperbolic parts of the critical manifold we know that there is a slow flow on a slow manifold 
that is 0(e)-close. Next, we recall an analog of the result for periodic orbits of the fast 
subsystem. The idea is to find a flow that approximates the flow on the family of periodic 
orbits. Consider the fast system 

^=x' = f(x,y) (16) 

such that f[T5"j) has a continuous family of periodic orbits j y (t) for each value of y in some 
neighbourhood D Q of y = with period T y that is uniformly bounded so that there are 
constants T a , T b > such that T a < T y < T b . For simplicity we shall also assume that 
each orbit 7^(t) is asymptotically stable with respect to the fast variables. It seems plausible 
that the full fast-slow system should have solutions (x(t), y(r)) such that the fast motion is 
approximated by the family of rapid oscillating periodic orbits: 



x( T ) ~ ly (^) • 
Formally plugging this result into the slow equation yields 

y = g(iy( T / e ),y)- 

The idea is that the slow motion on the family of periodic orbits can be obtained by averaging 
out the fast oscillations. Hence we might consider 



V. r,(Y):=^- I g(':ylt).Y)<lt. (17) 



1 



It will be convenient to make a change of variable t = TyQ and to set Ty(#) = 7y(T(y)#). 
This transforms (TT7I) to 

Y= [ g(T Y (e,Y))d6. (18) 







Assume that the solution Y{r) with initial condition Y(0) = Y stays inside D for < r < 
7~i . Then a classical theorem shows that our averaging procedure really produces the correct 
result with an error of order O(e). 
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Theorem 5.2. (\2% Let x be sufficiently close to T Yo (0 ) for some 9 . Then there exists 
a function 6(r) that satisfies a differential equation of the form 

e6 = ± + 0(e). 

J-Y 

Furthermore the following estimates hold 

x(t) = T Y (9(r)) + 0(e), 
y(r) = Y(r)+0(e), 

for 0(e\ loge|) < r < T\. 

Therefore we observe that the averaged systems (IT7j) - (IT51) appear in the condition f[T5]) . 
This means that points on the periodic orbit contained in dN are slow exit or entry points 
if the averaged flow is transverse to dN; we should view this averaged flow as a "slow flow" 
on the family of fast periodic orbits. Hence we have provided an analog for the condition of 
slow exit and entry points on the critical manifold. 

The next generalization step is now obvious. Consider a general family of invariant sets 
(e.g. tori) in the fast susbsytem. We can again average over the invariant measure of this 
family in the case of periodic orbits; the transversality conditions of this averaged flow will 
be exactly analogous to the previous cases. In practical applications this scenario does not 
seem to be needed very often as it does require three or more fast dimensions or a family of 
fast subsystems with one additional free parameter beyond the y-variables to be generic. 



6 Examples 

6.1 The Van der Pol Equation 

We re-consider Example 13.31 The time reversed version of Van der Pol's [51 [6] equation is 

x' = x 3 /3-x-y, , lg - 
y' = ex. 

The critical manifold Cq of ffT9j) is given by 

C = {(x,y)eR 2 :y = x 3 /3-x}. 

Two fold points are located at p± = (±1, =1=2/3). They naturally split the critical manifold 
into three parts 

d = c n {x < -l}, c m = c n{-i < x < i}, c r = c n{x>i}. 

A singular periodic orbit for (1191) exists for e = consisting of concatenations of solutions of 
the fast and slow subsystems. We choose N as a compact annulus containing the singular 
periodic orbit as indicated in Figure [6j In this case, we have 

S = lnv(N, O ) = C n N. 
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Figure 6: Critical manifold Co (grey), singular periodic orbit (black) and singular isolating 
neighbourhood (dashed blue) for Van der Pol's equation (|T9l) . 

Proposition 6.1. N is a singular isolating neighbourhood and (N,dN) is a singular index 
pair. 

Proof. Note that dN R S consists of four points, two on C m and one point each on G\ and 
C r . Since (Al)-(A4) hold, Theorem (14.11) implies that all four points are slow exit points as 
the slow flow is transverse at each point and pointing outwards with respect to N. Therefore 
Theorem 14.11 implies that N is a singular isolating neighbourhood. Next, we can apply 
Theorem 13.101 to see that L = dN and so (N, L) is a singular index pair. □ 

A direct Conley index calculation for (N,L), in combination with the existence of a 
Poincare section, can now be used to show that for e > sufficiently the singular periodic 
orbit perturbs to a periodic orbit of the full system ( TT§|) ; see [24] for the detailed calculation. 

Remark: Proposition 16.11 is well-known in Conley index and serves as one of the basic 
examples how to apply the theory to show the existence of a non-trivial invariant set in a 
dynamical system. We emphasize here that with our characterization of slow exit points 
in Theorem I4.1[ the problem has been reduced to the minimum amount of work regarding 
the checking of theorems; our transversality condition of the slow flow is much easier to 
understand and check than the conditions of Conley's Theorem 13.71 
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6.2 A Bursting Model 



We consider a modified Morris-Lecar model first proposed by Rinzel and Ermentrout [28 
x[ = y-0.5(xi + 0.5) -2x 2 {x 1 + 0.7) 

4 = 1.15 (o.5 (i + tanh^j^)) 

2/' = e(k-xi), 

where k is a parameter and < e 1. We note that (1201 exhibits special periodic orbits that 
are examples of bursting oscillations; see [17] for more details. Terman [31] and Guckenheimer 
and Kuehn [13] investigated fl20|) further focusing on the deformation of the periodic orbits 
under parameter variation related to a phenomenon called "spike adding". We shall not 
discuss these results further but refer to the original references. The important point in the 
current context is that the periodic orbits play a key role in the dynamics. Our goal is to 
construct a singular isolating neighbourhood for (I20p . In Figure [7J we show a bifurcation 
diagram for the fast subsystem 

x[ = y - 0.5( Xl + 0.5) -2x 2 ( Xl + 0.7)- 0.5 (l + tanh (^j^)) (an-1), (9U 
x' 2 = 1.15 (0.5 (l+tanh(2y|i))-x 2 ) cosh {^) , ' 1 ] 

where we regard y as a parameter. The diagram shows the continuation of an equilibrium 
point which traces out a projection of the critical manifold Cq = {x[ = = x' 2 }. All 
continuation calculations have been carried out using MatCont [32] • The two fold points 
Pi } r are fold (or saddle-node) bifurcations of the fast subsystem and at these points normal 
hyperbolicity is lost. They are located at 

Vi « (-0.0337,-0.0207,0.1365)=: (x ljh x 2}l , yi ), 
Pr « (-0.2449,0.0832,0.0085) =: (xi )r ,x 2>r ,y r ). 

The normally hyperbolic parts of Cq are separated by the fold points into three branches 

C b = C n {y < y r }, C m = C n {y r < y < y t }, C u = C n{y> y t } 

representing the lower, middle and upper parts of the S-shaped curve in Figure [7J It is easy to 
check that Cb is attracting and C m is of saddle-type. At yu ~ 0.075658 we find a subcritical 
Hopf bifurcation of (|2~T1) . The upper part of the critical manifold C u is repelling for y < yn 
and attracting for y > yjj. The unstable periodic orbits generated in the Hopf bifurcation of 
(j2ip undergo further bifurcations as indicated in Figure The bursting periodic orbits for 
the full system contain the perturbation of a segment connecting p r to the stable periodic 
orbits in the fast subsystem at y — y r . 

We shall not discuss how the periodic orbits connect back to Cb (see [13]) but note that a 
compact set N that is a potential singular isolating neighbourhood can have fast subsystem 
periodic orbits in dN. We shall focus on 

N := {(x 1 ,x 2 , y) e R 3 : (xi, x 2 ) G K, -0.04 < y < 0.084} 



0.5 ( l + tanh 
- x 2 ) cosh 



xi + 0.01 
0.15 

xi - o.r 

0.29 



(xi - 1), 



(20) 
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Figure 7: Bifurcation diagram for the equilibria of f l2T|) . Note that the S-shaped curve of 
equilibria also represents the projection of the critical manifold into (x2, y)-space. 



0.09^ 




Figure 8: Continuation of periodic orbits generated in the Hopf bifurcation for y = yu- The 
smaller inner orbits (in green) are unstable and the larger outer orbits (in blue) are stable. 
Stability changes at a saddle-node of limit cycles (LPC, limit point of cycles, in red). The 
thick middle curve indicates part of C u . Note that only a discrete number of periodic orbits 
are shown from the numerical calculation. 
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for a suitably chosen compact rectangle K C M 2 so that iV contains p/ and p r . Note that 
this choice of N is only a first step in the analysis of bursting orbits and has to be refined 
to prove their existence; see the technical problems discussed in [HI [TO]. The key point is 
that dN n {y = 0.084} contains two periodic orbits 71^ for the fast subsystem. Figure M,^) 
shows the two periodic orbits. 
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Figure 9: (a) Two periodic orbits in the fast subsystem at y — 0.084. (b) Integrals given in 
(|22|) for the two periodic orbits depending on the parameter k. 

From Proposition 15 . 1 1 we deduce that points on jj for j = 1,2 are slow exit /entry points 
depending on the sign of the integral 

h := / 3 gMs), 0)ds = [ ' k - ( 7i (s)) lt fo (22) 
Jo Jo 

where Tj is the period of jj. Figure Mjo) shows the value of Ij for different values of 
k G [—0.3,-0.1] which are typical values within which bursting periodic orbits occur [13J. 
We find that in this parameter range the outer cycle 71 consists of slow entry points while 
the inner cycle 72 consists of slow exit points. Note that dN D {y = 0.084} also contains 
a point in Co at which the slow flow is transverse to dN. Furthermore we can choose the 
compact rectangle K in the definition of N so that there is only one more point in dN PI Co 
lying on Cb- At this point the slow flow is again transverse to dN. By Theorem 14.11 both 
points are slow exit or entrance points. Next, we apply Theorem 13.51 to conclude that iV is 
a singular isolating neighbourhood in this case. 

Remark: Note that it is much more complicated to construct a singular index pair due 
to the presence of the Hopf bifurcation point of the fast subsystem on Co- We postpone this 
question on how to modify iV to future work. 
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